Quantitative laboratory observations of internal wave reflection on ascending slopes 
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Internal waves propagate obliquely through a stratified fluid with an angle that is fixed with 
respect to gravity. Upon reflection on a sloping bed, striking phenomena are expected to occur close 
to the slope. We present here laboratory observations at moderately large Reynolds number. A 
particle image velocimetry (PIV) technique is used to provide time resolved velocity fields in large 
volumes. The generation of the second and third harmonic frequencies are clearly demonstrated in 
the impact zone. The mechanism for nonlinear wavelength selection is also discussed. Evanescent 
waves with frequency larger than the Brunt- Vaisala frequency are detected and experimental results 
agree very well with theoretical predictions. The amplitude of the different harmonics after reflection 
are also obtained. 
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PACS numbers: 47.55.Hd Stratified flows. 47.35.+i Hydrodynamic waves. 



I. INTRODUCTION 

The oblique propagation of internal waves follows from 
the dispersion relation that monochromatic perturba- 
tions of frequency uo have to satisfy 



uj = ±iVsin0, 
where N is the Brunt- Vaisala frequency 



N = 



g_dp 

Po dz J 



(1) 



(2) 



g being the gravity, p(z) the ambient density profile and 
po a reference density. This dispersion relation shows 
that for a fixed frequency, the direction in which energy 
propagates with respect to the horizontal, 0, is fixed. 
Moreover, Eq. determines that phase and energy 
propagate in perpendicular directions. For set-up with 
Brunt- Vaisala frequency N independent of z, observa- 
tions of internal waves have invariably showed this 
transverse and oblique propagation. 

The above dispersion relation is obtained, away from 
any turbulent portions of the domain, by substituting 
a plane wave Aexp[z/c(xsin^ + zcosO)] of wavenumber 
k and amplitude A in the wave equation governing the 
horizontal velocity u 



Uttx 



■ Uttzz -N u x 



(3) 



where subscripts denote partial derivative, x and z refer 
to Cartesian coordinates and t to time. The vertical ve- 
locity w, the perturbation pressure p, the streamfunction 
ip and the perturbation density satisfy the same hyper- 
bolic equation ©. 



The presence of a solid horizontal bottom, a free sur- 
face or an oblique slope results in a reflected wave with 
the same intrinsic wave frequency as the incident one. 
Linear internal wave reflection on a sloping bottom has 
been treated analytically by different authors 0, E3 • 
The striking consequence of the geometric focusing of 
linear internal waves has also been reported [HI . Il2l Il3q . 
Under appropriate conditions, it leads to internal wave 
attractors in confined stably stratified fluids. For criti- 
cally incident waves for which the slope and energy prop- 
agation angles coincide, the linear inviscid analysis be- 
comes singular and an infinite amplitude of the reflected 
wave was predicted [T^ . 

Experiments on internal wave reflection were first per- 
formed by Cacchione and Wunsch 1] using conductivity 
probe measurements. A tidal-like excitation was gen- 
erated in a 5 m long tank by a horizontally oscillating 
paddle; the incident and reflected waves from a 15° slope 
were separated using periodogram estimates to compute 
wave amplitude and wavenumber. Although these pio- 
neering results were of poor quality (compared to what 
is possible nowadays), their shadowgraph experiments 
showed nevertheless striking microstructures along the 
slope, reminiscent of an array of vortices along the slope's 
boundary layer. Thorpe and Haines j24) measured dye 
band displacements on a 20° slope providing qualitative 
agreement with linear theory and also noticed three di- 
mensional boundary layer structures. Subsequently, Ivey 
and Nokes Q estimated the mixing efficiency above a 30° 
slope submitted to modal excitation and visualized with 
rainbow Schlieren technique: the weakening of the back- 
ground stratification was measured and the correspond- 
ing change in potential energy compared to the mechani- 
cal work provided by the wave-maker. The internal wave 
reflection mechanism has also been studied in close con- 
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nection with its importance for the creation of nepholoid 
layers by McPhee and Kunze ^i). More recently, Daux- 
ois, Didier and Falcon Q performed Schlieren experi- 
ments on critical reflection, focusing on the boundary 
layer upwelling events and providing good qualitative 
agreement with the previous weakly nonlinear study Q. 
Finally, using synthetic Schlieren measurements, Peacock 
and Tabaei [jjj showed the second harmonic generation. 
To our knowledge, no local quantitative amplitude mea- 
surements of internal wave reflection have been so far 
performed. 

Recently, the weakly nonlinear theoretical issue has 
been put forward. It has been shown by Dauxois 
and Young that the singularity can be healed using 
matched asymptotic expansion. Their analysis describes 
the build up of the reflected wave along the slope for a 
incident plane wave. Always from the theoretical point 
of view, the reflection of internal waves was revisited by 
Tabaei, Akylas and Lamb [2^ for the case of a narrow 
incident beam but for non-critical angles. In this case, 
they have predicted the generation of harmonics in the 
steady regime; this result has not been addressed exper- 
imentally yet. Besides, fully nonlinear numerical simula- 
tions have also been performed to examine the behavior 
of large-amplitude internal gravity waves impinging on a 

slope II ED, El 111 El 

Here we present the results of laboratory experiments 
in which a beam of internal waves is reflected on an 
oblique slope. Experiments were carried out in the 13m 
diameter Coriolis platform, in Grenoble, filled with salted 
water. The large scale of the facility allows us to strongly 
reduce the viscous dissipation along wave propagation 
and, moreover, quantitative results are obtained thanks 
to high resolution Particle Image Velocimetry (PIV) mea- 
surements. 

The paper is organized as follows. In Sec.|n] we present 
the experimental set up and discuss all control parame- 
ters. In Sec. IIIII we show the experimental results. We 
explain in Sec. IIII Al the spectral analysis used to dis- 
tinguish the different harmonics. In the following sec- 
tion IIII Bl we discuss the mechanism of wavelength se- 
lection. Section ITU CI is devoted to the evanescent waves, 
while Sec. IIII Dl discuss amplitude measurements. Sec- 
tion concludes and gives some perspectives. 



II. EXPERIMENTS 
A. Experimental setup 

We developed an original internal wave exciter inspired 
by Ivey et al. 8] in order to produce a two and a half 
wavelengths beam. A PVC sheet was compressed on both 
sides by seven squared arms (Fig. [IJ fixed on a long cen- 
tral axis via eccentric wheels. The axis was set in rotation 
by a stepping motor, generating a longitudinal oscilla- 
tion motion of 8 cm amplitude along the 60 cm width of 
the paddle. The frequency of excitation, proportional to 



the angular speed of the motor could be precisely moni- 
tored in order to vary the angle of propagation of internal 
waves. The wavelength was also varied from 11.3 to 12.6 
cm. The paddle itself was inclined at an angle <£> = 13° 
with the horizontal to increase excitation efficiency. If the 
main part of the energy is indeed transmitted to a beam 
of internal wave of frequency u, higher harmonics nuo are 
also excited by the oscillating paddle. Nevertheless, as 
their frequencies are higher, their propagation angles are 
also larger. Taking advantage of this fact, a screen has 
been appropriately located above the bottom-end of the 
glass plate so that all harmonics are reflected to the left, 
and do not perturb the region of interest. The emitted 
plane wave hits a 2 x 3 meters glass plane (see Fig. [Q , 
back-painted in black to avoid parasite laser beam reflec- 
tions. 




FIG. 1: Picture of the experimental setup. The 3 meters long 
oscillating paddle generates incident internal waves which are 
impinging on the inclined glass plane visible on the left. One 
also clearly sees the screen which avoids the propagation to- 
ward the glass plane of harmonics also generated by the pad- 
dle. This picture was taken before the 13 m diameter tank 
has been filled with 1 m of stratified salt water. 

The above experimental set up was put in the 13 m 
diameter Coriolis tank filled from below with salt water 
and stratified by computer controlled volumetric pumps 
from two 75 m 3 tanks, one filled with salt water and the 
other one with pure water. A fast conductivity probe 
and a temperature probe were lowered slowly into the 
tank using a controlled vertical micro-stepping motor to 
measure the stratification. The density probe was cal- 
ibrated using an Anton Paar densitometer accurate to 
0.0001 kg/m 3 and 0.01°C. The linearity of the result- 
ing density gradient was of very good accuracy, and only 
the upper 5 cm and the bottom 10 cm were not linearly 
stratified. Observations have thus been performed in the 
intermediate region where the Brunt- Vaisala frequency 
is a very well defined constant. 

The stratification of 2.2% over the 1 meter depth of 
the tank led to a Brunt- Vaisala period 2ttN~ 1 = 13.6s. 
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FIG. 2: Sketch of the vertical cross-section of the experimen- 
tal setup. The sinusoidal thick line on the top left indicates 
the position of the oscillating paddle, while the grey region de- 
fines the incident beam impinging on the sloping glass plane. 
The optic fibre, indicated by the solid line on the left, shines a 
laser light which is reflected on the underwater 45 degree mir- 
ror, making a thin vertical laser sheet necessary for the PIV 
measurements. Finally, the dashed line shows the position 
of the screen which avoids the propagation toward the glass 
plate of harmonics also generated by the paddle. A slit in the 
screen allows the laser sheet to pass through. All distances 
are in centimeter. 



All experiments we discussed were performed without ro- 
tation of the tank, and thus involve pure internal waves. 

We used the Particle Image Velocimetry (PIV) facility 
of the Coriolis Platform to obtain top and side views of 
the velocity field. The fluid was seeded with 400 microns 
diameter particles polystyrene beads that were carefully 
prepared by a process of cooking which decreases slightly 
the density and successive density separations. One thus 
obtains a fiat distribution of densities matching that of 
the salt stratification. This process ensures that there are 
equal number densities of particles at each depth. A sur- 
factant was added to prevent the polystyrene beads from 
clustering. It has been shown that the relaxation time 
for particles to attain velocity equilibrium 16] is about 
0.02s which is much shorter than the characteristic time 
scale of the flow. As usual, the particles are considered 
as passive tracers of the fluid motion. 

With this method, motion is vizualized by illuminating 
particles with a laser sheet, which are followed by a dig- 
ital camera. The laser is intentionally kept out of focus 
(approximately 1 cm sheet width), enabling tracking of 
particles despite some cross sheet displacement. Veloc- 
ity fields within the plane of the laser sheet are obtained 
by comparing patterns in two subsequent image frames 
(taken 1 s apart). 

The 6 Watt green Argon laser was placed above the free 
surface, but the light was reflected on a 45 degree mirror 
placed in the water (see Fig.EJ). Thanks to a second un- 
derwater 45 degree mirror, the images were acquired by a 
1024*1024 pixels CCD camera also located above the free 
surface. The cross-correlation PIV algorithm designed 
by Fincham and Delerce [5( was used to convert the im- 
ages into vector fields (stored in the standard file format 
NetCDF). This algorithm provides good anti-aliasing and 



peak-locking rejection procedures. Our resolution went 
to subpixels distortions, corresponding to submillimetric 
displacements of the fluid. Typical maximal displace- 
ment in an image pair corresponds to 5 pixels, with a 
measurement precision of 0.2 pixel for an individual field 
(4% relative precision). This is mostly a random error, 
so the precision on averaged fields is higher. Finally, the 
analysis of the .nc files was performed with Mat lab soft- 
ware. 

Typical experimental runs lasted about 20 minutes. 
The first 10 periods were considered as an initial tran- 
sient. Data were thus only recorded throughout the sec- 
ond stage during which the steady regime was attained. 

B. Characteristics of the incident wave beam 

The set of chosen excitation periods (see Table |l} leads 
through the dispersion relation to different angles of 
propagation ranging from 13 to 25°. As the slope angle 
a is set to 22°, these cases allow us to analyze subcritical 
(0 < a), supercritical (0 > a) or critical (0 ~ a) reflec- 
tions. Finally, as the wavemaker is slightly tilted (13°) 
from the horizontal, one has to take into account that 
the wavelength slightly varied from one run to another, 
around a typical value of Aq = 12 cm. 



Run 


1 


2 


3 


4 


5 


T(s) 


60 


49 


41.5 


36 


32 


(deg.) 


13 


16 


19 


22 


25 


a (deg.) 


22 


22 


22 


22 


22 


A (cm) 


11.3 


11.7 


12.0 


12.2 


12.6 


Umax (cm/s) 


0.412 


0.566 


0.775 


1.000 


1.225 



TABLE I: Summary of experimental runs with all control 
parameters. T is the excitation period, the angle of energy 
propagation, a the angle of the slope defined in Fig. and A 
the wavelength of the incident beam. Umax is the maximum 
horizontal velocity that has been measured. 

Previous experimental studies of internal wave reflec- 
tion have reported the importance of viscous dissipa- 
tion 0, IE E| to explain the observed steady-state so- 
lution. Here, by contrast, the large scale of the experi- 
ment allowed to work at a large Reynolds number. By 
taking the wavelength and the velocity amplitude of the 
incident internal wave beam, one gets for the Reynolds 
number Re c± 100, which suggests that viscous dissipa- 
tion might be negligible during the propagation toward 
the slope of the beam itself. This is what can be verified 
in Figs. BI where one has used a synchronous detection-like 
idea. We have filtered the PIV signal at the excitation 
frequency u (see Sec. IIII K\ for additional details). As 
the energy is expected to obliquely propagate with an 
angle 6 with respect to the horizontal, the relevant quan- 
tities are the velocity fields v s (s, a) = u cos — w sin and 
v a (s, a) = u sin 0+w cos 0. Both have been measured and 
the first one is reproduced with false color in panel (a). 



The longitudinal section presented in panel (b) with a 
solid line emphasizes that the amplitude of the incident 
beam is constant. The dashed line which shows the veloc- 
ity field V(j along the longitudinal section of panel (a) is 
as expected vanishingly small. This is an important nec- 
essary condition to discuss quantitatively the reflection 
process. Moreover, the picture emphasizes that incident 
phase planes, which correspond to the same color, are 
parallel to the direction of the incident group velocity, 
the latter being indicated by the dashed line. This is 
a clear demonstration of the orthogonality of the group 
velocity and the wave vector. Finally, panel (c) presents 
v s along the solid line of panel (a) and reveals that the 
width of the beam contains two wavelengths. This is an 
important point to draw a comparison with theoretical 
predictions derived for plane waves |3 . Again, the dashed 
line attests that the velocity field v a is vanishingly small. 




FIG. 3: Definition of variables. Three different coordinates 
systems will be successively used: the gravity-oriented one 
(x,z), the slope-oriented one (x\ z') and finally the incident 
beam oriented one (s, a). The PVC sheet compressed period- 
ically along the dotted line provides a shear flow orthogonal 
to this line, and thus tilted at a constant angle <I> with the 
vertical. This flow induces waves propagating at an angle 
with the horizontal, determined by the dispersion relation 
while a is the slope angle. 
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III. EXPERIMENTAL RESULTS 

A. Spectral analysis 

An important characteristic of the reflection of internal 
waves is the generation of different harmonics, theoret- 
ically predicted a long time ago |1 |2jJ, but only very 
recently experimentally observed jig. However the am- 
plitudes of the different harmonics can be very different 
and, consequently, hardly distinguishable even though 
their propagation angles are different. A typical example 
is presented in Fig. [5Ji. 

This is the reason why we developed a Fourier temporal 
analysis of the results, with filtering at the fundamental 
and higher harmonics frequencies. Given the two compo- 
nents of the velocity field provided by the PIV analysis, 
we compute the filtered velocity fields. For the horizontal 



FIG. 4: Picture in color. False-color pattern of the along- 
beam v s (x,z) velocity for — 16° and a = 22° in a vertical 
(x,z) section. The shaded triangle corresponds to the re- 
gion below the glass plate, unperturbed by the internal waves. 
The maximum velocity in this figure is 0.07 cm.s -1 . Panels 

(b) and (c) show the along-beam and cross-beam sections re- 
spectively, indicated by the dashed and solid lines in panel 
(a). The longitudinal variation v s is shown in panel (b) and 

(c) with solid lines, while the dashed lines correspond to the 
transversal variation v a . The longitudinal coordinate s and 
the transversal one a are defined in Fig. 

one u(x,y,t), one thus defines the different quantities 
2 f fl 

(u) n = / u(x, z, t) cos(nu;t + (/)) dt. (4) 

tl ~ ^ J to 

In the time interval [to, ti], the amplitude of the n-th har- 
monic (u) n is of course a function of the spatial variables 
x and z, but also of a constant phase 0, chosen with re- 
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FIG. 5: All panels in color. False-color velocity pattern in 
the vertical (x,z)-plane for the critical run 4 (see Table GJ. 
Panel (a) presents the instantaneous horizontal velocity field 
u(x,z,t) while panel (b) shows the phase- averaged velocity 
(u)i. Panels (c) and (d) show respectively the second (it;) 2 
and third (it;) 3 harmonics of the vertical velocity, (w) n , in the 
case of run 4 (see Table ^J. The two white arrows define the 
impinging region of the incident beam. In panels (c) and (d), 
the rays at the left of this region should thus not be taken 
into account: they have been generated by the screen. The 
maximum velocity in panels (a) and (b) is 2 mm.s -1 and, in 
panels (c) and (d), 0.4 mm.s -1 . 



spect to the fiat position of the sinusoidal paddle. This 
procedure is equivalent to band pass filtering the PIV 
time series at each point in the domain. 

Panels (b), (c) and (d) of Fig. show that this pro- 
cedure is an excellent tool to distinguish the different 
harmonics. Indeed, for this almost critical reflection case 



(run 4, a = 22°, = 22 ± 1°), it is clearly apparent in 
panel (b), that the reflected beam at frequency uo is ab- 
sent except a small slightly supercritical alongslope ray. 
Nevertheless, panels (c) and (d) show the emitted sec- 
ond and third harmonics propagating with steeper an- 
gles. The angles of propagation with respect to the hor- 
izontal for the second harmonic beam is O2 = 48 ± 1° in 
agreement with the theoretical value 

n = sin -1 (n sin (5) 

for n — 2. The third harmonic is evanescent; its char- 
acteristics will be discussed in Sec. 1111 CI Both harmon- 
ics are generated in the finite domain region, adjacent 
to the bottom, where the incident beam hits the slope. 
One should not take into account the few rays located at 
the left of the arrows (see Figs. 03 and 01): they were 
generated by the oscillations of the screen. The emitted 
second harmonic between the two arrows can still be spa- 
tially distinguished from these artifacts. It is important 
to emphasize that the color scales differ by a factor 5 
between the two first panels and the two last ones. 

In summary, Figs. exemplified that, even though the 
second and third harmonics are almost invisible from the 
instantaneous velocity field, they are very clearly appar- 
ent after the filtering procedure. This filtering method 
is consequently appropriate even when the amplitude is 
very small. To our knowledge, even if it was previously 
predicted 0, [2^ this is the first case where the third 
harmonic has been observed, and even more important, 
this is the first report of experimental quantitative mea- 
surements. 

It is important to notice that the vertical velocity field 
is presented in panels (c) and (d), rather than the hori- 
zontal one as in the first two panels. Indeed, as the sec- 
ond and third harmonics propagation angles are much 
steeper, the horizontal velocity field is of lower quality. 

An alternative possibility to keep the spatial confor- 
mation of the reflection process while distinguishing the 
different harmonics is to consider the specific kinetic en- 
ergy density field of each harmonic which can be deduced 
from Eq. Q as 

{E) n (x,y,<f>) = ±[{u) 2 n + (w) 2 n \. (6) 

An example is shown in Fig. for the subcritical run 2. 
One clearly distinguishes the incident beam impinging on 
the slope, and reflected downslope. It is clear that such 
energy plots give less contrasted results than the velocity 
ones as in Figs. [5] They are nevertheless extremely useful 
for amplitude estimates along cross sections of the beam 
as discussed below. 



B. Mechanism of wavelength selection 

The spectral decomposition presented in the previous 
section allows precise wavelength measurements across 
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FIG. 6: Color online. False-color energy pattern. The specific 
kinetic energy density of the first harmonic, {E)i, is shown in 
the vertical (x, z)-plane for the subcritical run 2 (see Table OJ. 
The maximum value is 5.10 -3 cm 2 .s -2 . 



the different beams, even in the impact zone. As the 
reflection surface is expected to play a key role in wave- 
length selection via boundary effects, we have focused 
our study in the boundary region along the slope. 

In Fig. Ufa), the contourplots of the filtered first har- 
monic (incident) and of the third harmonic (emitted) are 
superimposed in order to show the alongslope wavelength 
selection. It is clearly visible that the distance between 
the emitted phase lines has been strongly reduced com- 
pared to the incident ones. We emphasize that the wave- 
length selection mechanism appears to occur along the 
slope, where the superposition of the incident and the 
reflected beams generates nonlinear interactions. Assum- 
ing that this inner region plays a key role in the reflection 
of internal waves |3( , the only relevant dynamical behav- 
ior of the wave field has to be taken at z' = where the 
incident and the emitted waves can both be written as 
ijj(z' = 0) = Aexp[i(k x fX f — wt)], the amplitude A being 
different for the incident and the reflected waves. Nonlin- 
ear interactions may lead to second or third harmonics 
terms such as ^{z 1 = 0) = A2 exp[2i(k x 'x' — cut)] or 
^3(2/ = 0) = Aex.-p[3i(k x >x f — cot)]. The alongslope wave- 
length of the third harmonic is thus reduced by a factor 
three, as highlighted in Fig. |7fa) by the phase lines in- 
tersections with the slope. 

To gain further insight, Figs. [H| present quantitative 
measurements of this effect in the supercritical case 5. 
The cross sections of the first three harmonics at a fixed 
distance from the slope are presented. The difference in 
^-location are of pure geometrical origin and decrease, of 
course, when the cross-section is taken closer to the slope. 
These pictures are analogous to the theoretical results 
presented by Tabaei, Akylas and Lamb in Figs. 9(a-d) of 
Ref. H2. 

Let us stress that the above quantitative comparisons 
are possible despite the large difference between the in- 
cident and emitted amplitudes. The automatic elimina- 
tion of wrong vectors in the measurements by the high 
quality cross-correlation PIV algorithm (designed by Fin- 
cham and Delerce 15]) is the key point here to achieve this 
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FIG. 7: Pictures in color. Superposition of the vertical ve- 
locities contourplots of the incident first harmonic (normal) 
and of the emitted third harmonic (bold). Panel (a) presents 
the critical case (run 4) while panel (b) shows the subcritical 
case (run 2). The solid lines show the phase lines of the inci- 
dent first harmonic (w) 1 while the dashed ones correspond to 
the phase lines of the emitted third harmonic (1^)3- The two 
arrows in panel (a) indicate precisely the alongslope cross- 
section used in Fig. |H| 



goal. The measurements are indeed sufficiently precise 
to provide meaningful results even when the relative en- 
ergy ratio between successive harmonics is approximately 
0.1. This last ratio is of course directly linked to the pa- 
rameter £, used in Refs. 0, El f° r the small amplitude 
asymptotic expansion to describe theoretically the reflec- 
tion process. This is consequently not a limiting factor. 

The nonlinear wavelength selection can be clearly high- 
lighted using above the three pictures. Indeed, Fig. 
presents their spatial Fourier transform and clearly em- 
phasizes that the wavelength of the n-th harmonic is nA, 
explaining the wavevector tripling visible in Fig. |7fa) as 
theoretically predicted by Tabaei, Akylas and Lamb |2lj . 
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FIG. 8: Alongslope sections of the energy of the incident first 
harmonic (panel a), second harmonic (panel b) and of the 
emitted third harmonic (panel c) in the supercritical case (run 
5). The sections were obtained at 7 cm from the slope, as 
indicated by the arrows in Fig. ITT a). Note the difference in 
units for the ordinates. 



Experimental results in subcritical cases have unex- 
pectedly revealed an apparent different mechanism for 
the wavelength selection, when the slope angle a is larger 
than the angle of energy propagation. A typical ex- 
ample is shown in Fig. [Tfb). The second harmonic has 
disappeared and the wavelength of the third harmonic is 
equal to the incident wave when it is projected along 
the slope. This possibility was already mentioned by 
Thorpe [2^ since third order nonlinear interaction along 
the slope (z' = 0) may lead to third harmonics of the 
form ^s(z f = 0) = Ae-x.-p[i(k x >x' — 3o;t)], where the tem- 
poral frequency is tripled while the spatial frequency is 
kept constant. This is to our knowledge the first ex- 
perimental evidence of this possible nonlinear interac- 
tion. The transition from supercritical (third harmon- 
ics alongslope wavevector tripling and second harmonics 
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FIG. 9: Spatial Fourier transform of the alongslope sections 
presented in Figs. The solid line corresponds to < E >i 
(Fig. 0a)), the dashed one to < E > 2 (Fig. 0b)) and the 
dotted one to < E >3 (Fig. 0c)). This picture emphasizes 
the frequency tripling. 



alongslope wavevector doubling) to subcritical case (third 
harmonics alongslope wavevector conservation and ab- 
sence of second harmonics), is visible in Fig. Ufa) where 
for the critical case, both wavelengths are present at the 
frequency 3a; (same alongslope wavelength on the left of 
the impact zone, tripled on the right). However, this 
transition remains unexplained. 



C. Evanescent harmonics 

The angles of propagation of the different harmonics 
have been measured on the filtered patterns. All the re- 
sults are listed in Table [H] for the five runs of Table [I] 
The values are also plotted in Fig. ^| as a function of the 
pulsation u and of the harmonic number n. Panel (a) 
attests that the first and second harmonic are in perfect 
agreement with what is theoretically expected: all cor- 
responding symbols are on the straight line u = NsinO. 
This is not the case for the third harmonic in the three 
last runs. When the sine of the propagation angle is plot- 
ted versus the number of the harmonic as in panel (b), 
one also sees that the third harmonic is not aligned with 
the first two ones for runs 3 to 5. 



Run 


1 


2 


3 


4 


5 


0i (deg.) 


12 


16.5 


19.5 


22.5 


25 


2 (deg.) 


26 


35 


41 


48 


54 


3 (deg.) 


43 


52 


?65 


85 


81.5 


sin#i 


0.21 


0.28 


0.33 


0.38 


0.42 


sin #2 


0.44 


0.57 


0.66 


0.74 


0.81 


sin #3 


0.68 


0.79 


0.91 


1.00 


0.99 



TABLE II: Experimentally measured angles of propagation 
for the three different harmonics. The precision of the mea- 
sures is smaller than half a degree. The characteristics of the 
five runs are given in Table |H 




n 



FIG. 10: Angular dispersion for the different harmonics. 
Panel (a) presents the sine of the different angles of prop- 
agations versus the wave frequency uo. The different symbols 
corresponds to the first (triangles), second (plus) and third 
harmonic (stars). The solid line corresponds to the theoret- 
ical law uj = NsinO for 2ttN~ 1 = 13.6s. Panel (b) presents 
the same data as a function of the number of the harmonics 
for runs 1 to 5, respectively plus, stars, circles, triangles and 
squares). Error bars are smaller than the symbol size. 



As discussed for example by Tabaei et al |2l[ , the max- 
imum incident angle for which the n-th harmonic can 
propagate is sin _1 (l/n). For larger angles, the corre- 
sponding harmonic will be evanescent since the frequency 
would be larger than the Brunt- Vaisala top band fre- 
quency. For > sin -1 (1/3) ^19°, the third harmonic is 
thus found to be evanescent : this is the case for runs 3, 
4 and 5 of Table [H] explaining the three symbols not on 
the line in Fig. irUT a). In these cases, the third harmonic 
generated in the impact zone cannot propagate, and is 
thus trapped along the slope. This is what is shown in 
Fig. m Evanescent modes were also found experimen- 
tally 22] in a different context, nonlinear non-resonant 
interaction between two internal wave rays. 

If the evanescence of the wave is clarified, an angle of 
propagation, not vertical, can be measured as shown for 
run 5 in Fig.^J Another case is also visible in Fig.[5^d). 
It is possible to theoretically explain this angle as follows. 

Considering the linear equation © valid for inter- 
nal waves within the Boussinesq approximation, let us 
look for streamfunction solutions -0, evanescent in the z'- 
direction, i. e. orthogonally to the slope, but propagating 
in the x and z direction (see Fig.|3|for the definitions of 
these variables). Introducing the two components of the 
wave vector (k x ,k z ) and <5, the attenuation (or evanes- 



20 40 60 80 

x(cm) 

FIG. 11: Color online. False-color velocity pattern. The 
vertical third harmonic < w >3 is shown in the (x, z)-plane 
for the run 5. The maximum velocity is 0.4 mm.s -1 . 

cence) length, we thus look for solutions as 

z, tHfo e^ x + k * z ~ <"*) e~ z '/ 5 (7) 
=ipo e (ikx^in & / 5) %+(ik z — cos a / 5) z—icut ^ 

since z' = z cos a — x sin a. Introducing the above ansatz 
in Eq. Q leads to a complex equation. Separating real 
and imaginary parts and defining 7 = oj/N, we get 

o cos 2 a 7 2 — sin 2 a , N 

k z = /c^tana (l - 7" 2 ) . (10) 

For 7 > 1 , we can thus define the angle of propagation of 
the evanescent wave ev = atan(& x /A^). Equation (|T0j) 
leads to 

/1 / \ f cotan a\ . . 

ev (a,7) = atan I - _ 1 . (11) 

Figure[EI attests that the experimental results agree very 
well with formula ((11)) for run 4 and run 5. Unexpect- 
edly the case = 19° for which the third harmonic is at 
the limit of evanescence is however not described by the 
above model. Note also that an alternative theoretical 
description is proposed in Ref. [2lj . 

D. Amplitude measurements 

The determination of wave amplitudes is complicated 
by the beam-like appearance of the displacement field. 
Indeed, local temporal spectral analysis is not sufficient 
to provide reliable energy measurements, as phase and 
group speed are orthogonal. Spatial integration across 
the beam are needed to accurately evaluate the amounts 
of energy involved in each harmonics. 

First to validate the method of measurement, the mean 
displacement amplitudes A = \]v? + w 2 /uj of the inci- 
dent beam is plotted as a function of the incidence angle. 
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the difference between the second and third harmonic is 
small. 



FIG. 12: Angle of propagation of the evanescent third har- 
monic for run 4 and run 5. 



The wavemaker being tilted at a constant angle = 13° 
with the horizontal, it induces a shear motion perpendic- 
ular to itself (see Fig. EJ). The propagating part of this 
motion is thus expected to be proportional to sin(# + 
The very good agreement shown in Fig. justifies a 
posteriori the method. 




-5 
- a (deg) 



FIG. 14: Velocity amplitudes of the first three harmonics as 
functions of the difference between the emission angle and 
the slope angle a, for the five runs of Table |l] Stars, triangles 
and crosses correspond to the first, second and third harmonic 
respectively. 




sin(6+<I>) 

FIG. 13: Experimental measure of the displacement ampli- 
tude of the emitted wave as a function of the emission angle 
for the five runs of Tabled 

Determining the importance of the different harmon- 
ics emitted after the reflection process is an important 
issue to understand and hence describe theoretically the 
reflection process. After exploring several possibilities, it 
appears that the best method to quantitatively charac- 
terize the measurements is the following one. First, it is 
important to distinguish the two components of velocities 
vectors, parallel to the slope and orthogonal to it. The 
second one vanishes clearly close to the slope, satisfying 
thus the expected boundary conditions. Away from the 
slope, because of the possibilities of evanescence it is very 
delicate to get reliable amplitudes. On the contrary, the 
along slope component of the velocities is easier to deal 
with. It has no reason to vanish close to the slope and 
definitely did not. We have thus measured the ampli- 
tude of the wave by determining the maximum value of 
the different harmonics. Results are collected and shown 
in Fig. [21 The velocity amplitudes for the second and 
third harmonics are 4 times smaller than for the first 
harmonic, even for run 4 and run 5 when the third har- 
monic is evanescent. It is important to emphasize that 



IV. CONCLUSION 

In this paper, we have reported quantitative laboratory 
measurements of the propagation of internal waves. Our 
experimental set-up produces an incident beam of high 
quality. Its large scale allows to reach large Reynolds 
numbers, such that the effects of dissipation on the prop- 
agation is negligible. This is not anymore true for the 
critical reflection mechanism, which involves a strong re- 
duction of the wavelength, hence increased viscous ef- 
fects. Thus, one obtains a steady regime compatible with 
the hypothesis of previous theoretical models 0, E[ • 

The wavevector of the frequency n-th harmonics, pro- 
jected along the slope, is found to be proportional to n 
in the supercritical case. This is in agreement with the 
theory |2l|. A different selection mechanism is however 
observed in the sub-critical case, for which the wavenum- 
ber is equal to the incident one. This is in contradiction 
with [2l[, but in agreement with a more simple analysis 
previously proposed by Thorpe [23| . 

Harmonics with frequency higher than N cannot prop- 
agate and remain trapped near the slope. We document 
their existence for the first time and explain their spatial 
structure. 

The long two-dimensional internal waves exciter pro- 
duces very weak transversal variations apparently jus- 
tifying two-dimensional predictions in the (x, z)-plane. 
Considering a rotating tank would require to address the 
much tougher, fully three dimensional equation and it is 
not a priori clear what will happen. Work along this line 
is in progress. 
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